Large-scale instabilities in a non-rotating turbulent convection 
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A theoretical approach proposed in Phys. Rev. E 66, 066305 (2002) is developed further to 
investigate formation of large-scale coherent structures in a non-rotating turbulent convection via 
excitation of a large-scale instability. In particular, the convective-wind instability that causes 
formation of large-scale coherent motions in the form of cells, can be excited in a shear-free regime. 
It was shown that the redistribution of the turbulent heat flux due to non- uniform large-scale motions 
plays a crucial role in the formation of the coherent large-scale structures in the turbulent convection. 
The modification of the turbulent heat flux results in strong reduction of the critical Rayleigh number 
(based on the eddy viscosity and turbulent temperature diffusivity) required for the excitation 
of the convective-wind instability. The large-scale convective-shear instability that results in the 
formation of the large-scale coherent motions in the form of rolls stretched along imposed large-scale 
velocity, can be excited in the sheared turbulent convection. This instability causes the generation of 
convective-shear waves propagating perpendicular to the convective rolls. The mean-field equations 
which describe the convective-wind and convective-shear instabilities, are solved numerically. We 
determine the key parameters that affect formation of the large-scale coherent structures in the 
turbulent convection. In particular, the degree of thermal anisotropy and the lateral background heat 
flux strongly modify the growth rates of the large-scale convective-shear instability, the frequencies 
of the generated convective-shear waves and change the thresholds required for the excitation of 
the large-scale instabilities. This study elucidates the origins of the large-scale circulations and 
rolls in the atmospheric convective boundary layers and the meso-granular structures in the solar 
convection. 

PACS numbers: 47.27.Te, 47.27.Nz 



I. INTRODUCTION 



Large-scale coherent structures in a non-rotating tur- 
bulent convection at very large Rayleigh numbers are 
observed in the atmospheric convective boundary lay- 
ers 0, S H H, H |E Q, & H EE IH El, in numer- 
ous laboratory experiments in the Raylcigh-Benard ap- 

paratus M M EE III 13 EEJMJE1 M m 

and in direct numerical simulations |2J, |25|. Spatial 
scales of the large-scale coherent structures in a turbu- 
lent convection are much larger than turbulent scales 
and their life-time is larger than the largest time scales 
of turbulence. In the atmospheric shear-free convec- 
tion, the structures (cloud cells) represent large, three- 
dimensional, long-lived Benard-type cells composed of 
narrow uprising plumes and wide downdraughts. They 
usually embrace the entire convective boundary layer (of 
the order of 1-3 km in height) and include pronounced 
convergence flow patterns close to the surface. In the 
sheared convective flows, the structures represent large- 
scale rolls (cloud streets) stretched along the mean wind 
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Coherent structures in convective turbulent flows were 
comprehensively studied theoretically experiment ally 
and in numerical simulations LJ- L3 - LJ - 1^ l a , l a btl S tlS l QOl 

El El El El El EE El El liiniEIilsSlMlll 

|2E 123,123. However, some aspects related to the origin 
of large-scale coherent structures in non-rotating turbu- 
lent convection are not completely understood. Hartlep 
et al. (2003) noted that there are two points of view 
on the origin of large-scale circulation in turbulent con- 
vection |24j| . According to one point of view, the rolls 
which develop at low Rayleigh numbers near the onset 
of convection continually increase their size as Rayleigh 
number is increased and continue to exist in an average 
sense at even the highest Rayleigh numbers reached in 
the experiments [2I. Another hypothesis holds that the 
large-scale circulation is a genuine high Rayleigh number 
effect El- 

Recently, a new mean-field theory of non-rotating tur- 
bulent convection has been developed |3(J, |3l| . This the- 
ory predicts the convective-wind instability in the shear- 
free turbulent convection that results in the formation of 
large-scale motions in the form of cells. In the sheared 
convection, the large-scale instability causes generation 
of convective-shear waves. The dominant coherent struc- 
tures in this case are rolls. It was demonstrated [H |^| 
that a redistribution of the turbulent heat flux due to 
non-uniform large-scale motions plays a crucial role in 
the formation of the large-scale coherent structures in 
turbulent convection. 

In this study a theoretical approach |30l l3ll ] is devel- 
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oped further to investigate the formation of the coherent 
structures in the non-rotating turbulent convection. In 
particular, we investigated how the modification of the 
turbulent heat flux due to non-uniform large-scale mo- 
tions affects the critical Rayleigh number (based on the 
eddy viscosity and turbulent thermal conductivity) re- 
quired for the excitation of the convective-wind instabil- 
ity. We performed numerical study of the convective- 
wind and convective-shear instabilities in order to deter- 
mine key parameters that affect formation of the large- 
scale coherent structures in the turbulent convection. 

The paper is organized as follows. In Section II we 
discussed the physics of the formation of the large-scale 
coherent structures and formulated the mean-field equa- 
tions which describe the formation of the coherent struc- 
tures. In Section III we determined the critical Rayleigh 
number required for the excitation of the convective-wind 
instability in a shear-free turbulent convection. In Sec- 
tion IV we studied numerically the convective-shear in- 
stability in a sheared convection. Finally, conclusions are 
drawn in Section V. 



II. TURBULENT HEAT FLUX AND 
MEAN-FIELD EQUATIONS 

In this Section we discuss a redistribution of the turbu- 
lent heat flux due to the non-uniform large-scale motions 
as a key mechanism for the formation of the large-scale 
coherent structures in turbulent convection. Here we also 
formulate the mean-field equations which describe the 
formation of the coherent structures. Traditional the- 
oretical models of the boundary-layer turbulence, such 
as the Kolmogorov-type local closures, imply the follow- 
ing assumptions. Fluid flows are decomposed into or- 
ganized mean motions and turbulent flow. Turbulent 
fluxes are assumed to be proportional to the local mean 
gradients, whereas the proportionality coefficients (eddy 
viscosity, turbulent temperature diffusivity) are uniquely 
determined by local turbulent parameters. For example 
[32|. the turbulent heat flux reads F = (su) = -k t V5, 
where Krp IS the turbulent temperature diffusivity, S is 
the mean entropy, u and s are fluctuations of the veloc- 
ity and entropy, respectively. This turbulent heat flux F 
does not include the contribution from anisotropic veloc- 
ity fluctuations. 

Actually the mean velocity gradients can directly af- 
fect the turbulent heat flux. The reason is that ad- 
ditional essentially anisotropic velocity fluctuations are 
generated by tangling the mean-velocity gradients with 
the Kolmogorov-type turbulence due to the influence of 
the inertial forces during the life time of large turbulent 
eddies. The Kolmogorov-type turbulence supplies energy 
to the anisotropic velocity fluctuations which cause for- 
mation of coherent structures due to the excitation of a 
large-scale instability |30t l31j . Anisotropic velocity fluc- 
tuations are characterized by a steeper spectrum than 



the Kolmogorov-type turbulence 



steeper spectrum 1 

slum mm. 



The theoretical model 30] of the anisotropic velocity 
fluctuations and their effect on the turbulent heat flux 
includes the following steps in the derivations: applying 
the spectral closure, solving the equations for the second 
moments in the k space, and returning to the physical 
space to obtain formulas for the Reynolds stresses and 
the turbulent heat flux. The derivation are based on the 
Navier-Stokes equation and the entropy evolution equa- 
tion formulated in the Boussinesq approximation. This 
derivation |30l | yields the following expression for the tur- 
bulent heat flux F = (su): 
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where To is the correlation time of turbulent velocity cor- 
responding to the maximum scale of turbulent motions, 
W = V x U is the mean vorticity, U = + U z is the 
mean velocity with the horizontal Uj_ and vertical U z 
components, a is the degree of thermal anisotropy, and 



F* 
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(2) 



is the background turbulent heat flux that is the sum of 
the contribution due to the Kolmogorov-type turbulence 
(described by the first term in Eq. J5J)) and a contribution 
of the anisotropic turbulence caused by the shear of the 
imposed large-scale mean velocity (z) (the so-called 
counter-wind heat flux described by the second term in 
Eq. Q), 



(3) 



is a generalized anisotropic turbulent temperature diffu- 
sivity tensor. For turbulent convection b = (3/2) (2+7), 
7 is the ratio of specific heats (e.g., 7 = 7/5 for the air 
flow) and e is the vertical unit vector. The equation for 
the tensor Ky was derived in Appendix A in Ref. |30j us- 
ing the budget equations for the turbulent kinetic energy, 
fluctuations of the entropy and the turbulent heat flux. 
The anisotropic part of the tensor Kij (described by the 
second term in the square brackets of Eq. ©), is caused 
by a modification of the turbulent heat flux by the buoy- 
ancy effects. Note that for a laminar convection b is set 
to be zero and the coefficient of the turbulent tempera- 
ture diffusivity n T should be replaced by the coefficient 
of the molecular temperature diffusivity. The parameter 
a in Eq. Q is given by 



l + 4£ 



2/3 



- 1 



(4) 



where l± and l z are the horizontal and vertical scales 
in which the background turbulent heat flux F*(r) = 
(s(x)u z (x + r)) tends to zero. The parameter £ de- 
scribes the degree of thermal anisotropy. In particular, 
in isotropic case when l± = l z the parameter £ = and 
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a = 1. For l± <C Z z the parameter £ = — 1 and a = —9/2. 
The maximum value £ max of the parameter £ is given by 
£max = 2/3 for a = 3. The upper limit for the parameter 
£ arises because the function F*(r) has a global maxi- 
mum at r = 0. Depending on the parameter a the small- 
scale thermal structures in the background turbulent con- 
vection have the form of columns or pancakes (sometimes 
they are called as the small-scale thermal plumes). For 
a < 1 the small-scale thermal structures have the form of 
columns (l± < Z 2 ), and for a > 1 there exist the pancake 
thermal structures (l± > l z ) in the background turbulent 
convection (i.e., a turbulent convection with zero gradi- 
ents of the mean velocity). 

The terms in the square brackets in the right hand side 
of Eq. JIJ result from the anisotropic turbulence and de- 
pend on the "mean" (including coherent) velocity gradi- 
ents. These terms lead to the excitation of large-scale in- 
stability and formation of coherent structures. In Eq. Q 
the terms with zero divergence are omitted, because only 
divF contributes to the mean-field dynamics. Neglecting 
the anisotropic turbulence term in Eq. recovers the 
traditional equation for the turbulent heat flux. 

The physical meaning of Eq. is the following. The 
first term ex —tq a F* div Uj_ in square brackets in Eq. 
describes the redistribution of the vertical background 
turbulent heat flux F* by the perturbations of the con- 
vergent (or divergent) horizontal mean flows Uj_. This 
redistribution of the vertical turbulent heat flux occurs 
during the life-time of turbulent eddies. The second term 
tx r (a + 3/2) (WxF*) in square brackets in Eq. QJ de- 
termines the formation of the horizontal turbulent heat 
flux due to "rotation" of the vertical background turbu- 
lent heat flux F* by the perturbations of the horizontal 
mean vorticity Wj_. The last term cx To(W 2 xF*) in 
square brackets in Eq. describes the formation of the 
horizontal heat flux through the "rotation" of the hori- 
zontal background heat flux F^j_ (the "counter- wind" heat 
flux in Eq. ©) by the perturbations of the vertical com- 
ponent of the mean vorticity W 2 . These three effects are 
determined by the local inertial forces in inhomogeneous 
mean flows. A more detailed discussion of Eq. JIJ is given 
in Sections III-IV. 

The counter-wind turbulent heat flux (in the direction 
opposite to the mean wind) is well-known in the atmo- 
spheric physics and arises due to the following reason. 
In the sheared turbulent convection an ascending fluid 
element has larger temperature then that of surround- 
ing fluid and smaller horizontal fluid velocity, while a 
descending fluid element has smaller temperature and 
larger horizontal fluid velocity. This causes the back- 
ground turbulent heat flux F^ = —tq F* V 2 U±\z) in 
the direction opposite to the background horizontal mean 
sheared fluid velocity XJ^\z). 

We use the mean-field approach whereby the small- 
scale turbulent convection is parameterized. This is a 
reason why we do not use explicitly thermal plumes in 
the consideration. The main reason for the appearance of 



the large-scale coherent structures is related to the mod- 
ification of the heat flux by the non-uniform mean flows. 
The thermal plumes contribute to the modification of the 
turbulent heat flux. To a some extent, the redistribution 
of the turbulent heat flux can be interpreted as a redis- 
tribution of the thermal plumes. 

In order to study the formation of the large-scale co- 
herent structures in a small-scale non-rotating turbulent 
convection we used the mean-field Navier-Stokes equa- 
tion and the mean entropy evolution equation (with the 
turbulent heat flux QJ) formulated in the Boussinesq ap- 
proximation. These mean-field equations yield the fol- 
lowing linearized equations for the small perturbations 
from the equilibrium, U = U z 
and S = S - 
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g- t +U^V y -v T AjW = -aV x U , (6) 
| + Vy \ S = -(V • F) - (V,S<°>) U , (7) 



where v T is the eddy viscosity, Aj_ = A — V 2 , and 
(F* -e)[10aA_L - (8 a -3) A] U 
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K T (A + bV 2 z )S, (8) 



the turbulent temperature diffusivity. In order to 
derive Eq. ©, the pressure term was excluded by cal- 
culating the curl of the momentum equation. Equa- 
tions ©-0 allows us to study the linear stage of the 
large-scale instabilities. The variables U, W and S de- 
scribe the large-scale coherent structures. In Section III 
we study a shear-free convection with U^ ^ = 0, and 
in Section IV we investigate turbulent convection with 
the background (equilibrium) large-scale velocity shear 
U(°)(z) = <rze y and the background mean vorticity 
W(°) = VxU (0) = -ae x . 



III. SHEAR-FREE CONVECTION 

Let us consider a shear-free convection (U < -°- ) = 0). 
In the shear-free regime, the large-scale instability is 
related to the first term cx —To a F* div in square 
brackets in Eq. for the turbulent heat flux [33, El • 
When dU z /dz > 0, perturbations of the vertical velocity 
U z cause negative divergence of the horizontal velocity, 
divU^ < (provided that divU = 0). This strengthens 
the local vertical turbulent heat flux and causes increase 
of perturbations of the local mean entropy and buoyancy. 
The latter enhances perturbations of the local mean ver- 
tical velocity U z , and by this means the convective-wind 
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instability is excited. Similar reasoning is valid when 
dU z /dz < 0, whereas divUj_ > 0. Then negative pertur- 
bations of the vertical flux of entropy leads to a decrease 
of perturbations of the mean entropy and buoyancy, that 
enhances the downward flow and once again excites the 
convective-wind instability. Therefore, nonzero div Uj_ 
causes redistribution of the vertical turbulent heat flux 
and formation of regions with large values of this flux. 
These regions (where div Uj_ < 0) alternate with the low 
heat flux regions (where divUj_ > 0). This process re- 
sults in formation of the large-scale coherent structures. 

The role of the second term oc r (a + 3/2) (WxF*) in 
square brackets in Eq. 0} is to decrease the growth rate 
of the large-scale instability for a > —3/2. Indeed, the 
interaction of perturbations of the mean vorticity with 
the vertical background turbulent heat flux F* produces 
the horizontal turbulent heat flux. The latter decreases 
(increases) the mean entropy in the regions with upward 
(downward) local flows, thus diminishing the buoyancy 
forces and reducing the mean vertical velocity U z and 
the mean vorticity W. This mechanism dampens the 
convective-wind instability for a > —3/2. The above 
two competitive effects determine the growth rate of the 
convective-wind instability. A solution of Eqs. J5j and 
(|7|) in the shear-free convection regime yields the follow- 
ing expression for the growth rate 7 of long- wave pertur- 
bations: 
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1/2 



(9) 



where the parameter (3 — (Iq K)~ 2 3> 1, Zo is the max- 
imum scale of turbulent motions, 9 is the angle be- 
tween the vertical unit vector e and the wave vector 
K of small perturbations, g is acceleration of gravity. 
The analysis of the convective-wind instability was per- 
formed in Refs. H^El only for a small square of Brunt - 
Vaisala frequency. In particular, Eq. © was derived 
in Refs. El for the case \N 2 \ < gF*r K 2 , where 
N 2 = -(g • V)5(°) < is the square of Brunt- Vaisala 
frequency. 

In the present study we consider arbitrary values of 
the Brunt- Vaisala frequency, and we investigate the ef- 
fect of the modification of the turbulent heat flux (due 
to non-uniform large-scale motions) on the critical effec- 
tive Rayleigh number required for the excitation of the 
convective-wind instability. We also study here the effect 
of the anisotropy of turbulent thermal diffusivity (caused 
by the buoyancy) on the critical effective Rayleigh num- 
ber. To this end we rewrite Eqs. © and in a non- 
dimensional form 
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FIG. 1: Effective Rayleigh number versus the aspect ratio 
L z /L± of the perturbations for two free boundaries and dif- 
ferent values of the parameter n — (dashed line); for 
fj, = 0.7 (dotted line); /1 = 5 (dashed-dotted line). Here a = 1 
and b = 5.1. The classical Rayleigh solution for the laminar 
convection (b — 0) with the two free boundaries is shown by 
solid curve. 



where V = Pr T £ A U is dimensionless velocity, the length 
is measured in the units of the total vertical size L z of 
the system, the parameter b = 3 (2 + 7)/2 describes the 
anisotropy of turbulent thermal diffusivity caused by the 
buoyancy effect (see Eqs. ©), TZa — 6A 3 £Pr T is ef- 
fective Rayleigh number based on the turbulent viscos- 
ity, v T , and the turbulent temperature diffusivity, k t , 
A = L z /lo, Pr T = v T j ' K T is the turbulent Prandtl num- 
ber, the parameter fi is given by 



M = 



4 a, 
15" 
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ST is the mean temperature difference between bottom 
and upper boundaries of the turbulent convection, the 
parameter a, = 2 g tq F*/uq and To is the reference mean 
temperature. The last term oc /1 in the right hand side 
of Eq. Qll|l determines the modification of the turbulent 
heat flux due to the non-uniform large-scale motions, and 
the parameter 



Pr T 
lla~ 



1/3 



4gr F* 
\N 2 \L 2 



has the meaning of the normalized heat flux, 
|A 2 | = gST/[T L z ] and TZa = \N 2 \L\j[v T k t \. 



where 



A. Solution for two free boundaries 

Let us consider the solution of Eqs. lfTU|) and (|TT|l for 

two free boundaries, using the following boundary condi- 
tions 



V 



ViV = S = 



for 
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We seek for a solution of Eqs. I|10 JI -I|11 [) in the form 
V, S cx sin(7r n z) exp(7 t — i Ki • r) , 

where n is the integer number and K± is the horizon- 
tal component of the wave vector. The critical effective 
Rayleigh number (at 7 = 0) is determined by the equa- 
tion 



(Kl + (b+l) it 2 n 2 ) (Ki + tt 2 n 2 ) 2 = K 2 



lZa c 



-fi (Pr T TZa 2 c ) 1/3 [(2a + 3) K]_ - (8a - 3) tt 2 n 21 



where the critical effective Rayleigh number, lZa Cl is 
based on the turbulent viscosity and the turbulent tem- 
perature diffusivity. 

In the case of /i = (i.e., there is no modification of 
the turbulent heat flux due to the non-uniform large-scale 
motions), the critical effective Rayleigh number is given 
by 



TZa r — 



(Kl + tt 2 n 2 f(K\ + (b + 1) tt 2 n 2 ) 



K 2 ± 



(13) 



The minimum value of the critical effective Rayleigh 
number for the first mode (n = 1) for b = is lZa c w 
657.5. This is the classical Rayleigh solution for the 
laminar convection with two free boundaries. The criti- 
cal effective Rayleigh number increases with the increase 
of the anisotropy of turbulent temperature diffusivity 
(see Eqs. (01) described by the parameter b. Indeed, 
for b — 3.9 the critical effective Rayleigh number is 
TZa c rj 2247 and for b = 5.1 it is TZa c w 2722. 

The modification of the turbulent heat flux due to 
non-uniform large-scale motions strongly decreases the 
critical effective Rayleigh number. Indeed, Fig. ^ shows 
the effective Rayleigh number versus the aspect ratio 
L z /L± = K±/K z = tan# of the perturbations for differ- 
ent values of parameter ji. The increase of the parameter 
fi causes strong reduction of the critical Rayleigh number 
(see Table 0. 
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FIG. 2: Effective Rayleigh number versus the aspect ra- 
tio L z /Lj_ of the perturbations for two rigid boundaries and 
different values of the parameter fi: fi = (dashed line); 
(X = 0.7 (dotted line); /1 = 5 (dashed-dotted line). Here a = 1 
and b = 5.1. The classical Rayleigh solution for the laminar 
convection (b — 0) with the two rigid boundaries is shown by 
solid curve. 




FIG. 3: Effective Rayleigh number versus the aspect ratio 
Lz/L± of the perturbations for one rigid and one free bound- 
aries and different values of the parameter fi: /1 = (dashed 
line); \i = 0.7 (dotted line); /1 = 5 (dashed-dotted line). Here 
a — 1 and b = 5.1. The classical Rayleigh solution for the 
laminar convection (6 = 0) with the one rigid and one free 
boundaries is shown by solid curve. 



B. Solution for two rigid boundaries 



effective Rayleigh number is determined by the equation 



Now let us consider the solution of Eqs. (|10fl and 1111) 
for two rigid boundaries. In view of the symmetry of 
this problem with respect to two bounding planes it is 
convenient to translate the origin of z to be midway be- 
tween the two planes. Then fluid is confined between two 
planes z = ±1/2, and we seek for a solution of Eqs. I|1(J|) 
and l|llf) satisfying the following boundary conditions 



V = V z V = S = 



for 



1 

z = ±- 
2 



(14) 



We seek for the solution of Eqs. IjlUfl an d (|11|) in the 

form V, S oc exp(7 1 + q z — i Kj^ • r), where the critical 



[Ki-(i + b ) q 2 ](Ki- q y 
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-fi (Pr T Ka 2 c ) 1/3 [(2a + 3) K\ + (8a - 3) q 2 ] 



(15) 



The problem is symmetric with respect to the two bound- 
aries so the eigenf unctions fall into two distinct classes: 
the even mode with vertical velocity symmetry with re- 
spect to the mid plane and the odd mode with vertical 
velocit y a symmetry. Following the procedure described 
in [37J Ha, H3I , we adopted the even solution which has 
minimum critical effective Rayleigh number. Our nu- 
merical analysis showed that the anisotropy of turbulent 
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Boundaries 






Case 


two free 


one free and 


one rij 


pd two ri; 


gid 




L z /L ± 


lZa c 


L z /L± 


TZa c 


L z /L ± 


lZa c 


Laminar flow: 


0.707 


657.5 


0.854 


1101 


0.994 


1708 


Turbulent flow: 














n = 


0.891 


2722 


1.096 


4023 


1.280 


5547 


H = 0.7 


0.613 


1076 


0.697 


1328 


0.754 


4218 


fj. = 2.0 


0.578 


344 


0.645 


420 


0.688 


1743 


fi = 5.0 


0.568 


98 


0.628 


120 


0.662 


549 



TABLE I: Critical effective Rayleigh numbers for different 
types of the boundaries. Here for the turbulent flow a = 1 
and b = 5.1. The case of laminar convection is presented in 
Table only for comparison with the results obtained for the 
turbulent convection. 



temperature diffusivity described by the parameter b in- 
creases the critical effective Rayleigh number. In partic- 
ular, for b = the critical effective Rayleigh number is 
lZa c w 1707.8 This is classical Rayleigh solution for the 
laminar convection with the two rigid boundaries. For 
fi = and b = 3.9 the critical effective Rayleigh number 
is TZa c w 4683, and for b = 5.1 it is 7Za c ~ 5547. 

The effective Rayleigh number versus the aspect ratio 
L z /L± of the perturbations for two rigid boundaries is 
plotted in Fig. Increasing of the parameter /j, decreases 
both, the critical effective Rayleigh number and the as- 
pect ratio L z /L± of perturbations (see Table[IJ). If fi > 5, 
the behavior of the effective Rayleigh number drastically 
changes, e.g., there are two local minima for the effective 
Rayleigh number. 

C. Solution for one rigid and one free boundaries 

Solution for one rigid and one free boundaries can be 
obtained from the results for two rigid boundaries us- 
ing the odd mode. We use the domain from z = 
(the free boundary) to z — 1/2 (the rigid boundary). 



The anisotropy of turbulent temperature diffusivity de- 
scribed by the parameter b increases the critical effec- 
tive Rayleigh number. Indeed, for b — the critical 
effective Rayleigh number is lZa c p=> 1101 (the classi- 
cal Rayleigh solution for the laminar convection). For 
(i = and b = 3.9 the critical effective Rayleigh number 
is TZa c fa 3359, and for b = 5.1 it is TZa c ps 4023. The 
effective Rayleigh number versus the aspect ratio L z /L± 
of the perturbations for the one rigid and one free bound- 
aries is plotted in Fig. |3| Increasing of the parameter /i 
decreases the critical effective Rayleigh number and re- 
duces the aspect ratio L z /L± of the perturbations (see 
Table HI). 

Therefore, for these three types of boundaries the mod- 
ification of the turbulent heat flux due to the non- uniform 
large-scale motions strongly reduces the critical effective 
Rayleigh number (based on the eddy viscosity and turbu- 
lent temperature diffusivity) required for the excitation 
of the convective-wind instability. We summarized the 
final results for the above three types of the boundary 
conditions in Table [I] The case of laminar convection is 
presented in Table [I] only for comparison with the results 
obtained for the turbulent convection. 



IV. SHEARED TURBULENT CONVECTION 

In this Section we consider turbulent convection with 
a large-scale linear velocity shear U^°^(z) — aze y . In 
a sheared turbulent convection the mechanism of the 
convective-shear instability (33, is related to the last 
term oc To (W 2 xF*) in square brackets in Eq. The 
generation of the potential temperature perturbations by 
vorticity perturbations plays the key role in this mech- 
anism. Indeed, in two adjacent vortices with the oppo- 
site directions of the vertical vorticity W z , the turbu- 
lent fluxes of potential temperature are directed towards 
the boundary between the vortices. This increases per- 



turbations of the mean potential temperature and the 
buoyancy, and generates the upward flow between the 
vortices. These vertical flows excite vorticity perturba- 
tions, and the convective-shear instability mechanism is 
sustained. 

Let us consider an evolution of perturbations with zero 
y— derivatives of the fields U, W, S. We seek for a solu- 
tion of Eqs. ©-(0 in the form oc exp(7 1 - i K • r). The 
growth rate of the convective-shear instability of long- 
wave perturbations is given by [30l l3l| 

7 oc g F* tI K 2 ((3 A sin 2 6>) 2/3 , (16) 
where A = a tq is the shear parameter, the parameter 
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FIG. 4: Large-scale coherent rolls formed due to the 
convective-shear instability and aligned along the sheared 
mean velocity U'°'(«). The instability results in generation 
of the convective-shear waves which propagate perpendicular 
to the convective rolls. 
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FIG. 5: Range of parameters (L z /L±; L/lo) for which the 
convective-shear instability occurs, for different values of the 
shear parameter A=0; 0.05; 0.1; 0.2 . Here a = 1 . 

= (IqK)- 2 > 1 and K = ^/K 2 + K 2 . The convective- 
shear instability causes formation of large-scale coherent 
fluid motions in the form of rolls (see Fig.0J aligned along 
the imposed mean velocity U' ). The instability can also 
result in generation of the convective-shear waves with 
the frequency 

ftoc VZgF*T 2 K 2 {0\ sin 2 60 2/3 , (17) 

which implies the wave-number dependence, Q cx K 2 ^ 3 . 
The convective-shear waves propagate perpendicular to 
convective rolls (see Fig. QJ. The analysis of the 
convective-shear instability was performed in Refs. |3(i 
l3l| only for a small square of Brunt- Vaisala frequency 
and zero y— derivatives of the fields U, W, S. This cor- 
responds to the convective-shear instability for a very 
small component of the wave number along the imposed 
mean shear (i.e., uniform perturbations along the large- 
scale shear velocity). In this case the growth rate of the 
convective-shear instability is maximum. 

In the present study we consider arbitrary values of 
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FIG. 6: Range of parameters (L z /L±; L/lo) for which the 
convective-shear instability occurs, for different values of the 
angle tp between the horizontal wave vector and the x-axis: 
tp = 0° (solid line); tp = 18° (dashed line); ip — 30° (dotted 
line); tp — 90° (dashed-dotted line). Here a = 1 and A — 0.1 . 



the Brunt- Vaisala frequency and perform the numerical 
analysis of the convective-shear instability for nonzero 
y— derivatives of the fields U, W, S. We consider the 
eigenvalue problem with boundary conditions. We seek 
for a solution of Eqs. JSJ-Q in the form oc *&(z) exp(jt — 
iJ£± ■ r), where the eigenfunction $f(z) and the growth 
rate 7 of the convective-shear instability are determined 
by Eqs. ©-0- The system of the ordinary differential 
equations for the eigenvalue problem is solved numeri- 
cally with the following boundary conditions: U = U" = 
fjw = w = S = at z = 0, and U' = W' = S' = 
at z = 1, where /' — df/dz. We also take into account 
that for a turbulent convection, the turbulent Prandtl 
number can be estimated as Pr 1 w 4/(1 + Pr) w 2.34 
with Pr = 0.71 (for air flow). The latter estimate follows 
from the balance equations for the turbulent heat flux, 
the entropy fluctuations and the turbulent kinetic energy 

Let us consider the thermally isotropic (a = 1) turbu- 
lent convection. Figure [S] shows the range of parameters 
(L z /Lj_: L/Iq) f° r which the convective-shear instability 
occurs for different values of the shear parameter A. Here 

L = 1/ \J Lz 2 + LJ 2 and we assumed that a* = 1 . The 
case A = in Fig. [5] corresponds to the convective- wind 
instability (the shear-free turbulent convection). Inspec- 
tion of Fig.Elshows that the increase of shear is favorable 
for the excitation of the convective-shear instability. In 
Fig.|H|we plotted the range of parameters (L z /L±; L/lo) 
for the convective-shear instability for different values of 
the angle tp between the horizontal wave vector and 
the x-axis. Increasing the angle tp prevents from the ex- 
citation of the convective-shear instability (i.e., reduces 
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FIG. 7: Growth rates of the convective-shear instability and 
the frequencies of the generated convective-shear waves ver- 
sus: L z /L± and L/Iq. Corresponding dependencies on the 
parameters L/Iq are given for different L z /L± and vice versa. 
Here a — 1 and A = 0.1 . 



the range of parameters for which the instability occurs) . 
In Fig. we plotted the growth rates of the convective- 
shear instability and the frequencies of the generated 
convective-shear waves versus L z /L± and L/Iq. The 
curves in Fig. [7| have a point L* whereby the first deriva- 
tive d'j/dK has a singularity which is indicative of bifur- 
cation. The growth rate of the convective-shear instabil- 
ity for very small y— derivatives of the fields U, W, S, is 
determined by cubic algebraic equation [3(j • Below the 
bifurcation point, the cubic equation has three real roots 
(which corresponds to aperiodic instability without gen- 
eration of waves). Above the bifurcation point, the cubic 
equation has one real and two complex conjugate roots. 
In this case the convective-shear waves are generated. 
The source of energy for these waves is the turbulence 
energy. 

Now we perform the detailed numerical analysis of the 
convective-shear instability in order to determine the key 
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FIG. 8: Range of parameters (L z /L±; L/Iq) for which the 
convective-shear instability occurs, for A = 0.1 for different 
values of the degree of thermal anisotropy a: a — 1 (solid 
line); a = 0.9 (dashed line); a — 0.8 (dotted line); a — 0.7 
(dashed-dotted line). 
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FIG. 9: Growth rates of the convective-shear instability and 
the frequencies of the generated convective-shear waves for 
different values of the degree of thermal anisotropy a: a = 1 
(solid line); a = 0.9 (dashed line); a — 0.8 (dotted line); 
a — 0.7 (dashed-dotted line). Here A = 0.1 and L/Iq = 23 . 



parameters that affect this instability. First, we study 
the effect of the thermal anisotropy a on the convective- 
shear instability. Figure |H1 shows the range of parameters 
(L z /Lj_: L/Iq) for which the convective-shear instability 
occurs, for different values of the thermal anisotropy a. 
In Fig. El we plotted the growth rates of the convective- 
shear instability and the frequencies of the generated 
convective-shear waves for different values of a. The de- 




FIG. 10: Growth rates of the convective-shear instability 
versus the angle ip between the horizontal wave vector and the 
z-axis for different values of the degree of thermal anisotropy 
a: a = 1 (solid line); a = 0.9 (dashed line); a = 0.8 (dotted 
line); a = 0.7 (dashed-dotted line). Here A = 0.1, the values 
L z /L± and L/lo correspond to maximum growth rates of the 
instability. 
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FIG. 12: Growth rates of the the convective-shear instability 
and the frequencies of the generated convective-shear waves 
versus L z /L± for different directions ip of the lateral back- 
ground heat flux: ip = -45°; ip = 0°; ip = 45°. (a). The 
growth rates of the instability for L/lo ~ 10; (b). The growth 
rates of the instability for L/lo = 20; (c). The frequencies of 
the generated waves for L/lo = 10 (solid line) and L/lo = 20 
(dashed line). Here A = 0.1 . 
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FIG. 11: Range of parameters (L z /L±; L/lo) for which the 
convective-shear instability occurs, for different directions ip 
of the lateral background heat flux: ip = —45°; ip = 0°; 
ip = 45°. Here A = 0.1 . 

crease of the degree of thermal anisotropy a increases the 
threshold in the parameter L/Iq required for the excita- 
tion of the convective-shear instability. Figure ITUI shows 
the growth rates of the convective-shear instability ver- 
sus the angle tp between the horizontal wave vector 
and the x-axis for different values of a. Here the val- 
ues L z / Lj_ and L/Iq correspond to the maximum growth 
rates of the instability. For a > 0.7 the growth rate of 
the convective-shear instability attains the maximum for 
ip m > 0° . An increase of the degree of thermal anisotropy 
a increases the angle (p m . In the thermally isotropic 
(a = 1) turbulent convection the angle ip m = 18°, while 
for a = 0.8 (i.e., £ « 0.92), the angle ip m « 10°. Note 



that according to the atmospheric observations, the ob- 
served angle between the cloud streets and direction of 
the wind is of the order of 10° — 14°. The calculated 
angle cp m is in compliance with these observations. Note 
that the convective rolls are stretched in the horizontal 
plane in the direction perpendicular to Kj^ and the shear 
velocity is directed along the y-axis. Inspection of Figs.El 
and 1101 shows that decrease of the parameter a reduces 
the growth rates of the convective-shear instability. In 
Figs. l8llUl we considered the case a < 1 which is of inter- 
est in view of the atmospheric applications. 

Next, we study the effect of the lateral background heat 
flux (determined by the third term oc ((F*xe) • V) W 
in the right hand side of Eq. |JSJ), on the convective- 
shear instability. We introduce the angle ip between the 
horizontal component of the background turbulent 
heat flux and x axes, where the total background heat 
flux is F* = (F]_ cos ip,F]_ sini/j,F*). The angle ip is 
determined by the boundary conditions in the horizon- 
tal plane (e.g., by the temperature gradient in the hor- 
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FIG. 13: Growth rates of the the convective-shear instability 
and the frequencies of the generated convective-shear waves 
versus L/Iq for different directions ip of the lateral background 
heat flux: ip = -45°; ip = 0°; ip = 45°. (a). The growth 
rates of the instability for L z /L± — 0.5; (b). The growth rates 
of the instability for L z /L± — 2; (c). The frequencies of the 
generated waves for L z / ' L±_ = 0.5 (solid line) and L z /L± — 2 
(dashed line). Here A = 0.1 . 



izontal plane). Figure ITT1 shows the range of parame- 
ters (L z /Lj_; L/Iq) for which the convective-shear in- 
stability occurs, for different directions ip of the lateral 
background heat flux F^_. In Figs. H2fH?l we plotted the 
growth rates of the convective-shear instability and the 
frequencies of the generated convective-shear waves for 
this case, where F]_/F* = 0.5. Note that the background 
mean vorticity due to the imposed large-scale shear is 
W(°) = VxU(°) = -ae x . This is the reason why there 
is no symmetry with respect to the Y Z plane of the large- 
scale shear, i.e., the contributions to the convective-shear 
instability caused by the positive and negative angles ip 
of the lateral background heat flux are different. In par- 
ticular, the range of the convective-shear instability in 
the presence the lateral background heat flux with the 
positive angles ip is wider than that for the negative an- 
gles ip (see Fig. Illjl , On the other hand, even for the 
negative angles ip the range of the convective-shear in- 
stability is wider than that in the absence of the lateral 
background heat flux. Note also that in the presence of 
the lateral background heat flux with the positive angles 
ip, the convective-shear waves are not generated. This is 



reason why we plotted in Figs.ll2N13b the frequencies of 
the generated convective-shear waves only for ip < 0. 

Note that there are three groups of parameters in this 
study of the large-scale coherent structures formed in a 
turbulent convection: 

(i) the external parameters: the value of shear a and 
the background heat flux F* = (F± cos ip,F± smip,F*); 

(ii) the parameters which determine the background 
turbulent convection: the degree of thermal anisotropy 
a, the correlation time tq — Iq/uq and the parameter 
a* = 2 3 t F*/ul; 

(hi) the parameters related to the characteristics of the 
large-scale coherent structures: the aspect ratio of the 
structure Lj_/L z , the minimum size L of the structure, 
L x = Lj_ cos tp, the characteristic time of the formation 
of the large-scale coherent structures oc 7" 1 , and the fre- 
quency 51 of the generated convective-shear waves. 

The parameters related to the characteristics of the 
large-scale coherent structures are determined in this 
study. The external parameters and the parameters 
which determine the background turbulent convection at 
the present level of analysis are treated as free param- 
eters. The external parameters are determined by the 
boundary conditions. The degree of thermal anisotropy 
a can be determined by the budget equation for the two- 
point correlation function for the velocity-entropy fluc- 
tuations. This parameter has been recently measured in 
a laboratory experiment in turbulent convection in air- 
flow [2^ . In the range of the Rayleigh numbers 10 7 — 10 8 
(based on the kinematic viscosity and molecular diffusiv- 
ity) this parameter varies within the range from 0.5 to 
2. The parameter a* and the correlation time To = Iq/uq 
can be determined from the budget equations for the tur- 
bulent kinetic energy and vertical turbulent fluxes of mo- 
mentum and the entropy. The turbulent correlation time 
To and correlation length Ip are measured in laboratory 
convection (see, e.g., Ref. |23|). 



V. DISCUSSION 

In the present study we investigated formation of large- 
scale coherent structures in a non-rotating turbulent con- 
vection due to an excitation of large-scale instabilities. 
In the shear-free turbulent convection, the cell-like struc- 
tures are formed due to the convective-wind instability. 
The redistribution of the turbulent heat flux due to the 
non-uniform large-scale motions causes strong reduction 
of the critical effective Rayleigh number required for the 
excitation of the convective-wind instability. The effec- 
tive Rayleigh number is based on the eddy viscosity and 
turbulent thermal conductivity. We also found that the 
critical effective Rayleigh number increases with the in- 
crease of the anisotropy of turbulent temperature diffu- 
sivity caused by the buoyancy effects. 

In the sheared turbulent convection, the roll-like struc- 
tures stretched along the imposed large-scale sheared ve- 
locity are formed due to the large-scale convective-shear 
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instability. This instability produces the convective-shear 
waves propagating perpendicular to the convective rolls. 
We studied numerically the convective-shear instability 
and determined the key parameters that affect the for- 
mation of the large-scale coherent structures in the turbu- 
lent convection. In particular, we found that the degree 
of thermal anisotropy and the lateral background heat 
flux strongly modify the growth rates of the large-scale 
convective-shear instability, the frequencies of the gen- 
erated convective-shear waves and change the instability 
thresholds. 

The results described in this study are based on the 
linearized mean-field equations, and therefore, they can- 
not describe detail features of the turbulent convec- 
tion observed in the numerous laboratory experiments 
EllIlllllIllIMiJEIlSlHEl and in direct numer- 
ical simulations |24L f25j | . In particular, we made the fol- 
lowing assumptions about the turbulent convection. We 
considered a homogeneous, incompressible background 
turbulent convection (i.e., the turbulent convection with- 
out mean- velocity gradients). The nonuniform mean ve- 
locity affects the background turbulent convection, i.e., it 
causes generation of the additional strongly anisotropic 
velocity fluctuations by tangling of the mean- velocity gra- 
dients with the background turbulent convection. We as- 
sumed that the generated anisotropic fluctuations do not 
affect the background turbulent convection. This implies 
that we considered a one-way coupling due to a weak in- 
homogeneity of the large-scale velocity. Thus, we studied 
simple physical mechanisms to describe an initial stage of 
the formation of large-scale coherent structures in a non- 
rotating turbulent convection. The simple model consid- 
ered in our paper can only mimic the real flows associated 
with laboratory turbulent convection. Clearly, the com- 
prehensive theoretical and numerical studies are required 
for quantitative description of the laboratory turbulent 
convection. 

In spite of this very simple model it reproduces some 
properties of the semi-organized structures observed in 
the atmospheric turbulent flows |3l| . The semi-organized 
structures are observed in the form of rolls (cloud streets) 
or three-dimensional convective cells (cloud cells). The 
observed angle between the cloud streets and the mean 
horizontal wind of the sheared turbulent convection is 
about 10° — 14°, the lengths of the cloud streets vary from 
20 to 200 km, the widths from 2 to 10 km, and convective 
depths from 2 to 3 km. The ratio of the minimal size of 
the structure to the maximum scale of turbulent motions 
L/Iq = 10—100. The characteristic life time of rolls varies 
from 1 to 72 hours. Rolls may occur over water surface or 
land surfaces 0,0- Our study yield the following param- 
eters of the convective rolls: L/Iq = 10 — 100, the charac- 
teristic time of formation of the rolls ~ tq /j varies from 
1 to 3 hours. The life time of the convective rolls is de- 
termined by a nonlinear evolution of the convective-shear 
instability. The latter is a subject of a separate ongoing 
study. We have shown that the maximum growth rate 
of the convective-shear instability is attained when the 



angle between the cloud streets and the mean horizontal 
wind of the convective layer is about 10° — 17° in agree- 
ment with observations. We also found an excitation of 
the convective-shear waves propagating perpendicular to 
convective rolls. This finding is in agreement with ob- 
servations in the atmospheric convective boundary layer, 
whereby the waves propagating perpendicular to cloud 
streets have been detected ^lj. In addition, the motions 
in the convective rolls have a nonzero helicity in agree- 
ment with predictions made in Ref. |40| . 

There are two types of cloud cells in the atmospheric 
shear-free turbulent convection: open and closed. Open- 
cell circulation has downward motion and clear sky in the 
cell center, surrounded by cloud associated with upward 
motion. Closed cells have the opposite circulation 0. 
Both types of cells have diameters ranging from 10 to 40 
km, they occur in a convective layer with a depth of about 
1 to 3 km and the characteristic life time of cloud cells 
is about several hours. Our analysis shows that the min- 
imum threshold value of the effective Rayleigh number 
required for the excitation of the large-scale instability is 
attained at L±/L z ~ 2 (see Figs. 1-3, dotted and dashed- 
dotted curves), is in agreement with numerous observa- 
tions. The ratio of the minimum size of the structure to 
maximum scale of turbulent motions L/Iq = 5 — 15. The 
characteristic time of formation of the convective cells 
~ ro/7 varies from 1 to 3 hours. Therefore, the predic- 
tions of the developed theory are in an agreement with 
observations of the semi-organized structures in the at- 
mospheric convective boundary layer. The typical tem- 
poral and spatial scales of structures are always much 
larger then the turbulence scales. This justifies sepa- 
ration of scales which was assumed in the suggested in 
the theory. Note that the applicability of the mean-field 
equations for study of turbulent convection was discussed 
in Ref. BJ. 

In our study we consider non-rotating turbulent con- 
vection and apply our results to the atmospheric convec- 
tive boundary layers, where the shear is usually caused 
by wind. The rotation of the Earth usually affects the 
hight of the atmospheric convective boundary layer. The 
rotation can also affect the longitudinal spatial structure 
of the cloud streets. Note that in astrophysical applica- 
tions the shear (or differential rotation) can in general be 
a consequence of anisotropies in rotating systems. Our 
study can be also useful for understanding the origin of 
formation of the meso-granular structures in the solar 
convection (see [i^ 1 ). 
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